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In a paper by J. A. Lewis (B. S.T.J. , Jfi, 1961, pp. 1259-1280) general 
formulas for the electrical admittance of a piezoelectric resonator, having 
essentially one pair of electrodes, were derived in terms of motional param- 
eters associated with the normal modes of vibration of the device. The logical 
extension of this work to a resonator with N electrodes is presented here. 
Expressions are given for both the admittance and impedance matrices 
of the resonator. These matrices are expressed in terms of motional param- 
eters associated with, respectively, (i) the normal modes of vibration with 
all electrodes connected together, and (ii) all electrodes left open circuited. 
The electrical equivalent circuit for the 2-port characteristics of the N 
electrode resonator is given for two particular examples. 

I. INTRODUCTION 

General formulas for the electrical admittance of a piezoelectric 
resonator having essentially one pair of electrodes were derived by 
Lewis. 1 These formulas are consistent with those derived earlier for 
special cases such as long bars and large plates (see, for example, 
Mason 2 ). In Lewis' work the admittance function is expanded about 
its poles in an infinite series. The residue at one of these poles deter- 
mines the strength of the contribution of the normal mode, associated 
with the pole, to the overall vibrational behavior of the resonator 
when it is driven at a frequency close to the natural frequency of the 
mode. Surprisingly, the work of Lewis seems to have seen little ap- 
plication, as far as can be judged, except for that of Lloyd and Red- 
wood, 3 and Byrne, et al. 4 



* Most of the work described here is based on part of the author's doctoral 
thesis (University of London, 1966). 
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With the current interest in multi-electroded resonators, such as the 
monolithic crystal filter, 5 it is pertinent to consider the logical exten- 
sion of the work of Lewis to the case of an arbitrary resonator having 
N electrodes. A discussion of this problem has previously been pre- 
sented by the author, 6 and also by EerNisse and Holland. 7 

Included in Section II of this paper are the basic equations govern- 
ing the piezoelectric resonator, presented here for completeness. 

In Section III various integral relations are derived for use in Sec- 
tion IV where the properties of the admittance and impedance matrices 
are investigated. The electrical equivalent circuits for two particular 
2-port configurations of the N electrode resonator are also derived in 
Section IV, in order to illustrate the application of the admittance 
and impedance matrices. 

A brief list of the principal symbols used in the text is given below. 

1.1 List of Symbols 

p Mass per unit volume. 

p. Mass per unit area of an electrode. 

Ui Particle displacement vector. 

S kt Strain tensor. 

T kl Stress tensor. 

r< Traction (stress vector). 

<p Electric scalar potential. 

E t Electric field vector. 

D t Electric displacement vector. 

c E iikl Elastic stiffness tensor (measured at constant electric field). 

e nij Piezoelectric constant tensor. 

e£„ Dielectric constant tensor (measured at constant strain). 

n f Unit vector normal to, and outwards from surface of body. 

<J>„ Electric potential on the pth electrode. 

Qj, Total charge on the pth electrode. 

B Volume of the body. 

A TJnelectroded area of the body. 

A p Area of the pth electrode. 

w Angular frequency. 



\ — 2 



The tensor components above are referred to orthogonal Cartesian 
coordinate axes x { . The comma notation is used to indicate differentia- 
tion, e.g. Di,i = dDi/dXj , and the repeated index summation conven- 
tion is used, e.g., D,,, = D lA -\- D 2 , 2 + D 3 . 3 • 
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II. BASIC EQUATIONS OF A PIEZOELECTRIC RESONATOR 

The equations describing the steady vibrations of a piezoelectric 
body are listed below. 
The equations of motion: 

p\u, + r„., = 0. (1) 

The divergence equation of electrostatics (for an insulator) : 

D t .i = 0. (2) 

The piezoelectric constitutive relations: 

T« = cf ikl S kl - e nii E n , (3) 

D m = e mkI S kl + e s mn E n , (4) 

where 

S h i = \{Uk.i + 1*1,1), (5) 

and 

E n = -<f> 2n . (6) 

The symmetry relations 

E E E E /iy\ 

e„,y = e„n , (8) 

eL = eL . (9) 

2.1 Boundary Conditions 

The boundary conditions for the resonator shown in Fig. 1 will now 
be discussed. 

On the unelectroded portion of the surface A 

r, =0, on A, (10) 

DiUi = e Q (Ei) extn, = 0, an A, (11) 

that is, no surface tractions and zero external electric field exist normal 
to the surface. Note that (11) is an approximation which in practice is 
usually valid for materials with large values of ef m /e . The driving 
electrodes are assumed to be very thin metallic conductors with infinite 
conductivity. Potential <p p and charge Q p exist on the electrode area 
A„ . External electrical connections to the electrodes will not be specified 
at present. We pause to note, however, that <f> = at some point ex- 
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Fig. 1 — Arbitrary piezoelectric resonator with N electrodes. 

ternal to the resonator. Since we have neglected the effects of the 
external potential distribution, this "earth" point only has significance 
in connection with the topography of the external electrical circuit. 
The latter is assumed to interact only with the currents I v and potentials 
P on the electrodes of the resonator. The mechanical properties of 
the electrode are assumed here to be nonexistent except for a surface 
mass density p, . The surface of the resonator beneath an electrode is, 
therefore, assumed responsible only for exerting a force consistent 
with maintaining the acceleration of the electrode. The boundary 
conditions at the electrode can therefore be written as 



Ti = Pt^Ui , 

<p = constant, 



on A t 
on A. 



and either 



= $. 



or 



f D&idA, = -Q p . 



(12) 
(13a) 

(13b) 
(13c) 



Note that the choice between (13b) or (13c) as a primary condition 
is unimportant. 



III. PROPERTIES OF SOLUTIONS 

As is well known, the solution of the equations reviewed in Section 
II for a practical case is a formidable problem, and it is often neces- 
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sary to resort to some approximate method of solution. In this paper, 
we will not discuss the methods for solving the equations, but rather 
the nature of the solution assuming that it has been found. 

Equations (1) through (13) can be expressed in terms of w, and <£, 
from which it follows that 

p\iti + cf ikl u kii i + e nii <t>, nj = 0, (14) 

e nki u k , in - el0 inm = 0, (15) 

subject to the boundary conditions 

cfikiUkjUj + e kii <t>, k ni = 0, on A, (16) 

e ik iu k ,,iii - tfrtaJh = 0. on A, (17) 



with 

and either 



or 



and 



L 



<j> = constant on A p (18a) 

4 = *„ on A p (18b) 

D,n. dA p = -Q p , onA p (18c) 



cfikiUk.iUj + e kii (f}, k ni = p,\Ui , on A p . (19) 

We now note that (14) through (19) become homogeneous when 3> p = 
for all p. This latter condition represents one of the eigenvalue problems 
associated with Fig. 1, namely, that of mechanical vibrations possible 
when all electrodes are connected directly to the reference point. Other 
eigenvalue problems associated with Fig. 1 include those where some 
electrodes are open-circuited (Q p = 0) and the remainder are short- 
circuited (<f) r = 0) (i.e., connected to the reference point). 

3.1 Reciprocal Theorem 

Consider two solutions of (1) through (13) denoted, respectively, 
by (X', u' { , <£') and (X", u'/, <j>"). The two solutions could be, for example, 
those associated with two different sets of forcing parameters at dif- 
ferent frequencies. 

From (1) we have 

[ pXW dB + f u'iT'il, dB = 0, (20) 

J B J IS 
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and from (2) 

f tfD'tU dB = 0, (21) 

Jb 

where B is volume of the body exclusive of the electrodes. Using the 
divergence theorem, (20) may be written 

f pX'W dB - f ShT'i'i dB+ f u' t T' t fa, dA = 0, (22) 
Jb •> b Ja 

and (21) 

- f +'«!%' dB + [ t'D'/rii dA = 0. (23) 

J B J A 

Subtracting (23) from (22), and substituting for the surface condi- 
tions given by equations (10) through (13) we have 

f pX'u'M' dB+Y,\ P.X'W dA - £ & V Q' P ' 

J B p = l •» A v p=l 

= [ (T'/jSU - EM') dB. (24) 

J B 

Equation (24) is still valid when the primed and double-primed quan- 
tities are interchanged. Using this fact, we have 

(X" - X')| f puW dB+ JT f Pa uW dA - 2 {&M' ~ K'Q'v) 

L.J B p=l JAp J J)-1 

= f [(T'i'iS'a - EW) - (T'iiS'ii ~ E'/D'd] dB. (25) 

Jb 

The quantity on the right-hand side is zero by virtue of the con- 
stitutive equations (3) and (4). 
Equation (25) then becomes 



(X" - XOFGiW) = E WW - *£'&). (26) 

p=l 

where 

V(uW) - [ miM'dB+ E f p.u'M'dA. (27) 

Jb p-i Ja, 

Equation (26) is a special case of the reciprocal theorem given by 
Lewis 1 and discussed by Love 8 for the purely elastic case. 
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3.2 Orthogonality oj the Eigensolutions 

Consider two eigensolutions, (\ <n) , w,- n) , </> (n) ) and (X <m) , u\ m) , 4> (m) ) 
of the same homogeneous boundary problem. That is, $ ( v m) = if 
$p n) = and Q ( r m> = if Q ( r n) = 0. Thus, for two solutions of the same 
eigenset 

t (*r<r - *W) = 0, (28) 

and from (26) 

V(u\M mi ) = 0, X (n) ^ X (m) . (29) 

Thus, two solutions of the same eigenset satisfy the orthogonality 
condition (29). Also we have from (24), the Rayleigh quotient for the 
eigenvalue X ,n) 

2 J H(u\ ni , <f> M ) dB 

X< " > = w - = "Wj • (30) 

where 

H( Ui ,<f>) = l(r„S„ -EM. (31) 

3.3 Expansion in Terms oj Eigensolutions 

The solution to the inhomogeneous boundary value problem in- 
dicated by Fig. 1 can be expanded in terms of any of the sets of eigen- 
solutions. These expansions are verj r important when electrical behavior 
is of prime interest. We will show here how the forced vibrational 
solution (X, Ui , <£) may be expressed in terms of two of the possible 
expansions, namely: (i) the eigensolutions (X S(n> , wf (n) , <£ S(n) ) which 
correspond to the normal modes of vibration of the resonator with all 
its electrodes connected to the reference point, and (it) the eigen- 
solutions (X° (n) , u° <n) , <j> 0{n) ) for the normal modes with all the electrodes 
open circuited. 

For expansion (i) we set 

«i = u\ o) + 2 a ( V w (32) 

n-1 

and 

<t> = 4> (0) + Z« ( V (-> ; (33) 
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and for expansion (it) 

u . = „<•> + £ 6 w t,? w (34) 

and 

= * (o) + tt ( V w , (35) 

n=l 

where (wj o) , # (o) ) is the solution to the boundary value problem of 
(1) through (13), as X -» 0. 

Since we have not specified the means by which the electrodes are 
connected to the external electric circuit we will allow the parameters 
$ p and Q p to be of the general form 

*, = $< 0) exp fat), Q, = Ql 0) exp (*>0. (36) 

^•> = | $i«> | exp O'O, #" = I Q? I ex P (M- ( 37 ) 

Although it is immaterial how the charges Q p and potentials $ P are 
set up in relation to the external circuit, Q p and $„ are of course not 
independent. 

The coefficients o (n) in the first expansion can be found by noting 
that the equations of motion (1) and boundary conditions (12) require 

pXiii' 1 = P E (X S(n) - X)a w «? ( "\ in B (38) 

n = l 

and 

PM\" = P. E (X- S(n) - X)a ( V ( "\ on A p . (39) 

71 = 1 

On multiplying (38) and (39) by uf (m) and carrying out the indicated 
integrations and adding we have: 

f pkul'W* dB+it [ pM 0) u° (m) dA v 
= f) a M (\ SM - X)[ f pu* w u? lm) dB+it, f P.u! M u? (m) dA v ■ 

n =i LJb P=l J a, 

(40) 
We note from (27) that (40) may be written 

\V(vPvS lmr ) = Eo w (\ JW - X)7( W f"V (m) ). (41) 
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From the orthogonality condition (29), all terms on the right are zero 
except the term in o (m> , giving 

X V(u<W) 



X' s<m) - X V(uV m) uf im) ) 

By a similar argument we have for the coefficients in the second ex- 
pansion 

(X° (m) - X) V(uV m) uV m) ) 

Remembering the definition of u\ 0) , wf ( "° and u° (m) we have the follow- 
ing as a consequence of the reciprocal theorem (26) : 



and 



since 



x '"" F(»JV W ) = -Z Q l v el *! lm \ (45) 







and ( 


gJCl = 


• 0. 




Using (44) and (45) with (42) and (43) 








„ <m 


) 


p=l 


AgC, 








(X S(m 


, _ x)x s« m 


'7(ltf f - 


V c " 


°) 


b tm 


) 


-X 2 <r*r 

P =i 






(x O(m 


, _ x)x o (m 


'F( W ? ( " 


V ( " 


°) 



(46) 



(47) 

From (46), we see immediately that the contribution of the *S(w)th 
mode (eigensolution) in the expansion (32) and (33) is dominant when 
X — » X s<m> , if $p 0) is held constant with frequency. We also note that 
the amplitude of a (m) depends on the charge on the electrodes when 
the resonator is executing free vibrations corresponding to the S(m) th. 
mode (i.e., with all electrodes short-circuited). 

Equation (47) shows similarly that the contribution of the 0{m) th 
mode is dominant in the expansion of (34) and (35) when X — » X 0(m) 
if Qp 0) is held constant. Also the amplitude 6 (m) depends on the po- 
tentials on the electrodes when they are open circuit with the resonator 
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executing free vibrations corresponding to the 0(m) th mode. When 
applying the expansions (i) and (ii) it should be realized that assump- 
tions have been made concerning the completeness of the eigensets. 

IV. THE ELECTRICAL ADMITTANCE AND IMPEDANCE MATRICES 

4.1 General Considerations 

The admittance matrix y Pa for the A r -electrode piezoelectric resonator 
is defined by 

h - Z Vv& ■ (48) 

Similarly z VQ , the impedance matrix, is here defined by 

<*>,= £z PQ I Q + R, (49) 

a = l 

where 

h = jo>Q v (50) 

and R is a constant depending on the external circuit configuration. 
The relationships (48) and (49) are postulated on the basis that the 
equations of the resonator are linear and that their use is restricted to 
steady vibrations. 

We will now derive various properties of y PQ and z„ . First we note 
from (13c) and the divergence theorem that 

Z /, = ico f; Q p = -j. E [ D f n 4 dA p = -jco f D,,,- dB = 0. 

p = l p-l p=l J A, J B 

(51) 

Equation (51) is simply Kirchhoff's current law, for the conservation 
of charge. We note from (48) and (51) that 

£ t V»*. - 0, (52) 

p-l 0=1 

and since $ is arbitrary, the sum of each column of the y pq matrix is 
zero, i.e., 

£ 2/p„ = 0. (53) 

„=i 

We will now use the reciprocal theorem to show that both y PQ and z PQ 
are symmetric. Consider the solutions for two sets of potentials $ p 
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and <pp' having the same frequency. Then from (26) 

f; wa» - WW) = o. (54) 

Using (48) and (54) 

L 2 {Vv&W - &.*£'*£) = o, (55) 

p=l 0=1 

and (49) and (54) 

i"Zt, MW - zM'Q'p) +RZ(Qp- Q") - o. (56) 

p=l 0=1 p=l 

Since <!>', and & p ' are arbitrary in (55) we must have 

V P * = y«v • (57) 



In (56) Q' p and Q' v are arbitrary and ^ Q p = 0, so 



•PQ ~QP 



(58) 



As a consequence of (57) and (53) 

E Z/ P o = 0. (59) 

= 1 

It has been shown in this section that the impedance and admittance 
matrices of an iV-electrode piezoelectric resonator have properties 
similar in many ways to those of an ^-terminal passive electrical 
network." 

4.2 Expansions for the Electrical Parameters 

We may now make use of the eigensolution expansions of (32)-(33) 
and (34) -(35) to inquire into the admissible forms for y p<1 and z pa as 
functions of frequency. 

For expansion (a) 

Q, = Qi 0) + t a c-) « c-) (60) 

m = l 

#p = ^ 0) + Z « w «! w = *J°\ all < (m) being zero. (61) 

m = l 

Then using (46) and (60) 

q v = q'p o) + tt *£:T*\x ■ ( fi2 ) 

o=i m=i (A — a; 
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where 

QS (m) /-)S (m) 
5 — ¥1 (63) 

x S(m, y(wf (m, wf (m) ) 

For expansion (b) 

Q P = Qi m) + Z b (m) Q°< m) = <ff\ all Q° {m) being zero, (64) 

m-l 

and 

•t-«4 rt + Z6 ( " , *; ( - > . (65) 

m = l 

So using (47) and (65) 

$ = $«" - F E pq ^" - , (66) 

where 

F (»> = 5e — *s (67) 

x o(m) F( W r m, ^ (m) ) 

We now define the charge-potential relations for the solution of the 
static boundary value problem (X = 0) as follows: 

Ql 0) = t CM" (68) 

a-l 

and 

•4" = E F l *Q? + R. (69) 

0=1 

It is assumed, from now on, that the static parameters such as C' VQ 
are such that quadratic forms like C*,*,*, are positive definite. Proof 
of this depends on energetic considerations. 

We now put 

X 0(m) = coL and \ aim) - a&. , (70) 

and define 

CS = c»«- Y.cS\ (7D 

■a— 1 

We obtain from (62) and (71) the admittance matrix of (48) in the 
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form 



y M - wcff + E 



coke' 



. = 1 (Wflm — CO ). 

From (66) and (69) the impedance matrix of (49) is of the form 



1 



2 Fi 



z = _ Ipw _ y 

jW I m = l (0U m — CO ) 



(72) 



(73) 



Restricting our interest for the moment to an element y pq of the 
admittance matrix, we observe that the form of (72) is analogous in 
form to the admittance of the electrical network in Fig. 2. However, 
from (59), which is valid for all frequencies, we require 



but from (63) 



so 



Z <#' = o, 



Q(m) > q 



E C2° < o. 



(74) 
(75) 
(76) 



Therefore, several elements of C^ may be negative. 

The form of an element of the impedance matrix (72) suggests an 
analogy with the circuit of Fig. 3 but, in view of the proceeding dis- 
cussion, it is again probable that several of the elements F™ are negative. 
It should be noted, however, that F™ > 0. 

4.3 Driving Point Functions 

From (49) the driving point impedance Z£ , at the two terminals 
p-q, when all other terminals are left open-circuit is given by 





for ypq.C^Cp^CmHCJJ and o, m = o, Rm 



Fig. 2 — Electrical network for admittance representation. 
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C( C 2 C m 





FOR z pq ,C = i/FgJ.Cm^i/F^AND w m =w Am 
Fig. 3 — Electrical network for impedance representation. 



(77) 



In terms of the expansions for z pg given by (73), we have from (77) 

i r °° 2 /?w i 

P ° JW L «=1 W^m — • J 



where 



and 



»c«) = A*p gg ' (79) 

X 0< " ,, F(^ (m, wf (m, ) 



F w = f w _ 2F £ +F w i (80 ) 

So, clearly F (m) > and therefore, the analogue circuit of Fig. 3 is 
"physical" for Z p " . 

At first sight, it would appear that one could easily derive a driving 
point admittance for the p-q port when all other terminals are shorted. 
In fact, it must be found by appropriate manipulation of either y pq 
or z VQ and, in general, the expression includes many elements of either 
matrix. We can, however, define a driving point admittance for the 
p-terminal, when all other terminals are connected to the reference 
point, i.e., 

Yi = y„ . (81) 

The analogue electrical circuit for y„ , namely Fig. 2, is again physical. 

4.4 "Black Box" Matrices for a Two-Port 

Before calculating any "black box" transfer matrices it is convenient 
to define a transformed admittance matrix valid for the resonator and 
its external circuit. In Fig. 4, the terminals of the resonator are all 
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interconnected, there being a physical component with admittance 
y* 9 connected between terminals p and q. The currents flowing into 
the iV-terminal network and resonator as a whole are 



where 



and 



k = z; g.*. . 



z/ Pa = y„ - y f 



Vvp = 



E 



y* 



(82) 



(83) 



(84) 



t/ p , is denned by (48) and y£ = y * , as can be seen from Fig. 4. We 
may now form two-port networks. For the purposes of further dis- 
cussion, any connections made externally to the two-port will be as- 
sumed to be consistent with 

/.' = -I' p , l' T = -I'„ , V, = * p - <!>, , V 9 = # 8 - $ r . (85) 

4.5 Electrically Symmetric Two-Port Resonator 

Further discussion, with all y£ finite for the iV-terminal resonator, 
will not be continued. The reduction of a iV-terminal network to a 
2-port is discussed by Weinberg. 7 

4.5.1 Two-Port With N-2 Terminals Shorted 

We will now consider a simple case where all terminals except p 
and q are connected directly to s, and an admittance y* 9 is connected 




JL 

Fig. 4 — External electrical connections to the resonator. 
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between p and q as shown in Fig. 5. We will also assume that the con- 
struction of the resonator is such that it is electrically symmetric 
with respect to the ports. We have under these circumstances 



where 



and 



/-' = v > y _j_ v > y 

1 p Uvv ' V I Uvi ' Q 

■* « = Vvq ' v i Uvv ' q i 



y'va = Vpi - Vt 



(86) 
(87) 

(88) 

(89) 



We now consider the electrical lattice network of Fig. 6 as an analogue 
of transfer characteristics of (86) and (87). The analogue (Fig. 6) is 
physical if Y a and Y b are realizable with physical components. 
Now 



and 



Y a = Vw - y PQ , 



Y b = Vw + Vm • 



(90) 



(91) 



Using (72) and (63), Y a and Y b can be expressed in terms of the eigen- 
solution expansion as follows: 



;:| '■" ;; r 4t (x stm > - X)F( W f ( "V (m) ) J ' 

W .C)g;<-) _, ( y 



i<\_( 



y, = H (C> + cj?) + Z 



=t (k Sim) -X)FK (m V <m) ) 



f, im) Q! (m) ) ~\ 



(92) 



(93) 



Ip 


p 


y P E q 


q 


u 




ft 






ft 




fvp 


) 


1 


V q 










Ip' 


s 




r 


I q' 



Fig. 5 — Two-port system for N-2 electrodes short-circuited. 
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Fig. 6 — Electrical analogue of Fig. 5. 



Also since we have taken y pp = y Qa , then from (63), either 

.e, = _ Q s M or gci = q sm (94) 

It therefore follows that the S(m)th eigensolution may only contribute 
to one of Y a and Y b , depending on sign of Q p (m) /Q^ (m) . We also see 
that the electric circuit of Fig. 2 is a physical analogue for both Y a 
and Y b . 

4.5.2 Two-Port With (N-4) Terminals Open Circuit 

The symmetrical resonator with (iV-4) terminals open circuit is 
shown in Fig. 7. We now use the z PQ matrix of (49), (67), and (73) to 
derive the impedance matrix of this two-port, again subject to the 
restrictions of (85). We find that 



where 



and 



V = 7 J A- 7. J 

' p ~* PP P ' PQ Q 

I g = 7 pll l v -f" 7pp' Q , 

& pa PQ I *jr *"pr ^aq • 



From (82) and (67) 



4 

p) L 



Z„„ - — I G, 
p> 



- z 



XG< 



=! (A° (m) - X)- 



(95) 
(96) 

(97) 

(98) 

(99) 
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Ip 
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q 


la 




~~ Ut— 


» 
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4^ 


V q 


Ip 
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r 


X q 



where 



Fig. 7 — Two-port system for 2V-4 electrodes open-circuited. 



rr(") = pw i pw _ 97T ,0> 
<Jp P * pp ~ * •> -* p« i 

u " " x 0(m, 7( M ? (^ V m, ) ' 

X «->7(ii? f - ) fi? < - ) ) 



(100) 
(101) 

(102) 
(103) 



Ujij — ' pq I l ar l VT ' »fl • 

Also since the resonator has been taken to be symmetrical, i.e., 

Z vp = Z Qq , (104) 

it follows, from 67, that 

(*; cm> - tf'"') = ±« (m) - *? (m) ). (105) 

If we represent the transfer equations (95) and (96) in terms of the 
lattice analogue of Fig. 8, subject to the restrictions of (85), we have 

Z a = -7- (G v ° p — G V " Q ) 
ju L 

(106) 
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Fig. 8 — Electrical analogue of Fig. 7. 

We note by virtue of (105) that the 0(m) th eigensolution only 
contributes to one of Z a and Z b depending on the sign of ($" (m) — $? (m) ) 
/(#? <m) - $ r 0(m) )- It also follows that each of Z a and Z b have the circuit 
of Fig. 3 as a physical analogue. 

V. DISCUSSION AND CONCLUSIONS 

It has been shown how the electrical behavior of a piezoelectric 
resonator with N electrodes, represented by an admittance or im- 
pedance matrix, can be determined from the eigensolutions for free 
vibrations of the resonator. The results for the admittance obtained 
by Lewis 1 are contained here as the special case N = 2 in (72). The 
impedance of the two-electrode resonator is described by (73). This 
result was not given by Lewis 1 since he did not consider the alternative 
expansion of the open circuit eigensolutions. It could be argued that 
since impedance is simply the reciprocal of admittance, residues of 
one could be found from the other. This would, however, involve 
rather cumbersome calculations. Furthermore, if an approximate 
theory is used to generate the eigensolutions, the expansions may only 
be valid in a small frequency range, thus making the calculation of, 
say, the residues of the impedance from the admittance expansion 
impossible. 

Returning to the general case of the N- electrode resonator, the elec- 
trical behavior can be predicted in terms of the admittance or im- 
pedance matrices of (72) and (82). If external electrical components 
are to be connected between the N electrodes, and a two-port com- 
posite network is to be formed, its transfer characteristics can be 
deduced from either matrix. In the particular case of the symmetric 
resonator with N-2 of its electrodes connected together, the admittance 
matrix provides simple results, whereas the impedance matrix is use- 
ful for the case of A r -4 open-circuited electrodes. 

Before concluding, it might well be asked what reasons there are 
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for preferring a motional parameter representation of the electrical 
characteristics over the direct method of determining the impedance 
or admittance matrices from the solution to the inhomogeneous bound- 
ary value problem. These are essentially twofold. Firstly, if a two- 
port iV-electrode reasonator is to be designed to realize some transfer 
function or driving impedance, or so on, an appropriate synthesis 
procedure will usually automatically yield these motional parameters, 
leaving only the task of physically realizing a resonator with these 
parameters! Secondly, if the inhomogeneous boundary value prob- 
lem is being solved, the inevitable numerical calculations are least 
likely to be accurate in just those ranges which are of prime interest, 
namely, the poles and zeros of the admittance or impedance matrices. 
Furthermore, considerable computing time would be lost, compared 
with the motional parameter method, if some transfer characteristic 
of the derived two-port were to be obtained for a large number of 
frequencies in a narrow band. 
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